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Abstract 

We consider black p-brane solutions of the low energy string action, computing scalar pertur- 
bations. Using standard methods, we derive the wave equations obeyed by the perturbations and 
treat them analytically and numerically. We have found that tensorial perturbations obtained via a 
gauge-invariant formalism leads to the same results as scalar perturbations. No instability has been 
found. Asymptotically, these solutions typically reduce to a AdS( p+2 ) x S^ 8 ~ p ^ space, which, in the 
framework of Maldacena's conjecture, can be regarded as a gravitational dual to a conformal field 
theory defined in a (p + l)-dimensional flat space-time. The results presented open the possibility 
of a better understanding the AdS/CFT correspondence, as originally formulated in terms of the 
relation among brane structures and gauge theories. 

PACS numbers: 04.50. Gh, 04.70.Bw, 04.30.Nk 
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I. INTRODUCTION 



String theory and the subsequent idea of branes have been, in recent years, the almost 
standard theory describing the physics of quantum space-time, especially near the Big Bang 
or even before it Jj. The discovery of the relation between anti-de Sitter space physics and 
Conformal Field Theories on the boundary of that space, the so-called AdS/CFT correspon- 
dence {2,2] implied further interest in the structure of the string- membrane theory. 

The p-brane extended solutions are considered fundamental in the understanding of the 
non-perturbative string theory regime. They interpolate AdS p+ 2 x S d ~ p ~ 2 and d- dimensional 
Minkowski space-time [4]. This connection was important for the conjecture presented by 
Maldacena in 1997 j^j , which opened the way for the gravitation- field theory dualities. In this 
context, a better understanding of the perturbative dynamics of the p-brane solutions are 
relevant for the structural aspects of the AdS / CFT correspondence and its latter extensions. 
Such extensions can provide new hints about Yang-Mills theory with special interest in what 
concerns the difficult question of a Quark Gluon Plasma, see for instance [5 7 ] . Besides, 
in the framework of AdS/CFT correspondence it is possible to study the glueball mass 
spectrum analyzing the dynamics of a scalar field in the near horizon limit of the black 
p-brane solutions [81-11 01] . The poles of the retarded function of the simplest glueball state, 
generated by the operator O = Tr(F 2 ), are the quasinormal modes of the dual AdS black 
hole in the corresponding near horizon limit. 

A fundamental feature of the p-brane backgrounds is the possible existence of event hori- 
zons. In this sense, they may be viewed as generalizations of the usual four-dimensional 
black holes. Perturbations of black hole solutions are well known [ill . [12I and several numer- 
ical methods exist, being under full control to handle the information gathered from such 
perturbations 
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15|. 



We intend here to first define a perturbation of a p-brane solution using standard sepa- 
ration of variables and subsequently treat, analytically and numerically, the wave equation 
for the scalar perturbation. The employed methods are largely independent, aiming to a 
cross-check of the results. We also consider gauge-invariant gravitational perturbations. The 
results turn out to be exactly the same as scalar case. 



One very recent work complement our analysis presented here 16( . But although the pre 



sented work and 16[ are complementary and relevant in terms of the AdS/CFT correspon- 
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dence, they treat different geometries and focus on different issues. The results presented in 



this paper address directly the role of the brane structure (in the sense presented in 2). 



1) 



on the gravitation-field theory duality, specifically searching for possible instabilities. 

The paper is organized as follows. Sec. II provides reviews p-brane background considered 
in this work. In Sec. Ill, the perturbative dynamics is formulated and developed, followed 
by Sec. IV and V where the non extreme and extremal scenarios are specifically treated. In 
Sec. VII some final comments are presented. 

II. p-BRANE SOLUTIONS 

Solutions of ten dimensional Supergravity describing the so-called p-branes are well 



known. 



by lairzj: 



*et us consider the bosonic sector of type II Supergravity in ten dimensions, given 



^2 

p+2 



(1) 



Einstein's equations with N electric charges and p dimensions is obtained 
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s = csoW / d * m ^ r (R + 4 (V0)2) - wms 

where l s is the string length, g the determinant of the metric tensor g a b, R the Ricci scalar, 
4> the dilaton field and F p+ 2 the field strength of the potential A p+ \. 

The solution o: 
from the Ansatz 

ds 2 = ds 2 w _ p + e a Y, d V ld y l > (2) 

8=1 

where ds\ Q _ p is the line element with lorentzian signature in (10 — p) dimensions, a is a 
function of x, that is the bulk's radial coordinate, and the meaning of N as a charge arises 
from the Gauss Law. We can write a full solution as 

ds 2 = -A(x)dt 2 + B(x) [dr 2 + r 2 ^_ x ] + C(x)dx 2 + x 2 D(x)dtt 2 8 _ p , (3) 

where A(x) = (1 - (a/x) 7 - p ) (1 - (b/x) 7 - p y 1/2 , C(x) = (1 - (b/x) 7 - p ) ai (1 - (a/x) 7 ^)" 1 , 
B(x) = y/l - (b/xy-p, D(x) = (1 - (b/x) 7 -n a2 , with 01 = -| - §Eg and a 2 = \ - {g}- 
The mass per unit volume is M = ^ 7 _ P ^ Kl [(8 ~ p) fl7 ~ p ~ b 7 ~ p ], the electric charge = 
i [ab]^ 7 ~ p ^ 2 , Ki = (2Tc) 7 d p l p , k 2 = d p g s l 7 s ~ p } and g s is the string coupling, l p the Planck 
length in ten dimensions and d p = 2 5-p 7r( 5-p ' )//2 r((7 — p)/2). Absence of naked singularities 
implies 



N 



M> ,„^ 1P+1 ■ (4) 



Considering the non-extreme scenario, the maximal extension of the metric describes 
a black brane geometry, with an event horizon located at x — a. If p ^ 3, a curvature 
singularity is present at x = b, while if p = 3 we observe that, in addition to the outer 
horizon at x — a, there is also an inner horizon at x — b, with the singularity at x — 0. 
That behaviour is observed in the the Kretschmann scalar JC p (x) = R a b c dR abcd , where R a bcd 
are the components of Riemann tensor, as seen in the expression for the divergent term, 



(i-«rr^> 



where 5(p) = [(1 + p) + 2(5 — p)} if p is even, and 5{p) = ^(p— 3) — l)(p — 5) 
+ §0-3)0-5) if pis odd. 

For extremal p-branes the metric reads 

cfe 2 = E{x) [-dt 2 + dr 2 + r 2 dn 2 p ^] + F(x)dx 2 + x 2 G(x)dnj_ p , (6) 

where E{x) = y/l - (a/x) 7 -P, F(x) = (1 - (a/a;) 7 - p ) 71 , G{x) = (1 - (a/x) 7 -P) a \ 7l = 
«i — 1. 

In the extreme case, the curvature singularity is located at r = a and the metric does not 
have an extension if p ^ 3. We have a curvature singularity, but its structure depends on 
the value of p. If p = 6 the singularity is time-like, and the proper definition of a Cauchy 
problem is delicate. On the other hand, if p = 0, 1, 2, 5, the singularity (r = a) is null jsj, and 
therefore much milder. In spite of the absence of a event horizon, the manifold is globally 
hyperbolic, and the wave problem is well-posed. For the extreme case and p = 3, there an 
analytic continuation of the metric beyond r = a and we have again a black hole solution 
as pointed out in 



III. SCALAR AND GRAVITATIONAL PERTURB ATIVE DYNAMICS 

We initially consider a massless scalar field in the background of our 10-dimensional 
solution. We will show in the following that this scenario is more general. This perturbation 
is described by the Klein-Gordon equation 

A 10 $ = [Ap(r, 00,-x)) + A 10 _ p (t, x, A(g_ p) )] $ = , (7) 
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where the first term refers to the subspace dr 2 + r 2 dQ 2 _ 1 and the second to the bulk coordi- 
nates (t, x, A(8- p )). We denote the angular coordinates in dQ 2 ^ and dfll_ p respectively by 
by 9(p-i) and A 8 _ p . 

Such equation can be separated by the Ansatz Q(x A ) = 
^2im R i( r ) Y im{0i)Y J L<,^ , L(t,x)Y L g(X j ), where Y lm {0i) and Y Lq (Xj) are the well known 



spherical harmonics in (p — 1) and (8 — p) dimensions respectively 
differential equations 



d 



r ( P -i) j r 



,(p-i) 



i\dRi 



dr 



P 2 - 



l{l+p-2) 



Ri = o, 



resulting in the 



(8) 



d 2 V 



where u(x) 



A{x) 
B(x) 



P 2 + 



ill 

B(x) 
x 2 D{x)' 



'' +4tA,*l + u(i)*i = 0. 



A(x) 



(9) 



LiL + 7 — p) . Moreover, (3 is a constant arising from the 



brane {r, and bulk {t, x, \(8- P )} variables separation. Moreover, A x is a differential 

operator given by 

£ Ua(x)B(x)C(x)D(*-p)x 8 -v£ 



A. 



(10) 



^/Al^)Bl^)C{x^IW^x 8 -P 
The solution of equation (JSj) is i$(r) = Air 1 "^ 2 J 7 (/3r) + A 2 r l ~ p/2 Y 1 (f3r), with 7 = 
|a/p 2 — 4p + 4 + 4/(Z + p — 2), Ax and A 2 being constants, J 7 (/3r) and Ky(/3r) the Bessel 
functions. Finiteness at origin implies A 2 = and -Rz(r) = A\r x ~ p l 2 J 7 (/3r). Therefore, /3 
has a continuous spectrum of allowed values, and we notice in (J9j) that the its square acts 
as a mass for the Klein-Gordon field. Performing the same analysis for a time independent 
scalar field in the near horizon limit of the metric fl3]), the /3 2 parameter can be interpreted 
as the glueball mass. 

A "time independent approach" can be explored expanding the function with a 

Laplace-like transform [19j . Within this approach, we obtain the equation 



d? 
dr': 



Z L + [k 2 - V(x)] Z L = 



(11) 



where we defined the tortoise coordinate as dr*/dx 



y/C(x)/A( X ), * L (t,x) 



J e tult b(x)ZL(x)du with b(x) 
tential is given by the expression 
'A(x) 



C (8-P)/2B( x )p/4£)( :r )(8-p)/4 J 



k 2 = uj 2 — (3 2 and the effective po- 



V(x) 



B(x) 



P 2 + 4Kn L ( L + 7 -p)-ih 

x 2 V(x) b[x) 



h(x)b(x) — g(x)b(x) 



(12) 



6 



where the primes denotes differentiation with respect to x, h(x) = A(x)/C(x), and g(x) 



A(x) d 



111 



A{x)B{x)(D(x)xf-P 
C(x) 



C(x) dx 

We can also consider the problem of the linear perturbations using the gauge-invariant 
formalism proposed by Ishibashi et al [121 ] . In this formalism we expand the gravitational 
perturbations in terms of tensor harmonics Ely, and perturbations of Einstein equations 
are expressed as a group of equations for gauge invariant quantities. Such quantities are 
grouped in three types: tensor, vector and scalar. For the sake of simplicity, we only 
consider in the following the tensor sector of gravitational perturbations. The spacetime 
is considered as describing an m + n-dimensional manifold A4, which is locally written as 
the warped product g a pdz a dz 13 = g a b(y)dy a dy b + f(y) r y i jdx i dx : ' , where jij(y) is the metric of 
an n-dimensional maximally symmetric space of constant spatial curvature, and g a b{y) the 
metric of an arbitrary m- dimensional space time. 

Following reference 12] the following equation for the gauge-invariant quantity Ht can 
be obtained: 

UH T + ^-Dr ■ DH T - l{l+ J~ p) H T = 0, (13) 

where □ is the D'Alembert operator written on the metric g a b{y)- Introducing in the above 
equation the master variable $ = f^Hx we found the same result that we have already 
obtained from the scalar Klein-Gordon equation. 

At this point it is appropriate to make the following important observation: the spec- 
trum of quasinormal frequencies for the scalar field perturbations contains extra modes with 
respect to the tensor perturbations, because the modes for the last case only appears for 
multipole numbers equal or greater than 2. Thus, for the black p-brane, we need only to 
consider a test scalar field perturbation. Extracting the I > 2 terms for the obtained spec- 
trum of scalar quasinormal frequencies, we obtain the spectrum for the tensor gravitational 
perturbations. 



IV. NON EXTREME CASE 



The effective potentials derived above determine the perturbative dynamics. Of particular 
importance for this dynamics are the quasi normal modes. They are defined as solutions 
of the wave equations which satisfy the in-going and out-going boundary conditions. These 
modes are particularly relevant in the intermediate time behavior of the perturbation. 
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With arbitrary L, two different and independent numerical tools will be used in this work 
to calculate the quasinormal frequencies: a "frequency domain" approach based on a sixth 
order WKB techniq ue |20(|. and a "time domain" method based on a numerical characteristic 



integration scheme 



2lM23|. Both algorithms are well established. 





Figure 1: Effective potential for several values of (5. The p-brane parameters of are a = 2, b = 1, 
L = 1; and p = 3 (left), p = 6 (right). 



The WKB expressions are usually accurate and straitforward. But the approach is not 
generally applicable. For instance, in Fig. [1] the effective potencial is presented for a few 
values of (3 with p = 3, 6. We observe that the maximum of the effective potential decreases 
as (3 increases for a given p. For a sufficiently large value of (3 the potential becomes 
negative. This behavior appears explicitly for p = 6 with (3 = 1. Therefore, we cannot 
obtain the quasinormal frequencies for all values of p and (3 using the WKB formula. The 



instability for effective potentials that exhibit a negative gap is not excluded 24|, |25[. Direct 
time integration can be used for such scenarios. We have found no instabillities after an 
extensive exploration with (3 2 > 0. 

Within the "time-domain" approach, we have observed the usual picture in the pertur- 
bative dynamics. After the initial transient regime, the quasinormal mode phase follows as 
well as a late-time tail. The tail phase is strongly dependent on the value of the parameter 
p. For (3 = 0, we have a non-oscillatory power-law decay. But if (3 ^ 0, the tail is oscillatory, 
with a power-law envelope. Typical profiles are shown in Fig. El 

Given the potential, we use the sixth order WKB technique to obtain the quasinormal 
frequencies k. From the numerical data Z L (t,Xf ixec [), it is possible to estimate the funda- 
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Figure 2: Log-log graph of the absolute value of Zi(t,x fi xe d)- The quasinormal and tail phases 
are indicated. The p-brane parameters are p = 0, a = 2, b = 0.5, L = 1 and (3 = (top), (3 = 1 
(bottom). 

mental quasinormal frequency with reasonable accuracy. Some results from both methods 
are given in Tables H] HT1 and [TTT1 for j3 = 0. The concordance between them is good. However, 
notice that for p = 6 and L = 0, our result should be taken with reservation. Higher over- 
tones are not accessible by the "time-domain" technique. The corresponding WKB results 
are presented in Table IIIII 



The dependence of the frequencies u = \Jk 2 + f3 2 on (3 was also investigated. Both WKB 
and direct integration methods were employed, although the time evolution approach is not 
applicable for large (3, since in this regime the massive tail dominates from a very early 
time. Nevertheless, it should be reliable for small /3. Generally, we observed that for large 
values of the mass parameter, as /3 increases the frequencies becomes more oscillatory and 
less damped. One intriguing point was seen in a specific choice of parameters, namely a = 2, 
b = 0.5, L = and p = 2. In this case, the WKB and time evolution methods give discrepant 
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Table I: Fundamental quasinormal frequencies with a = 2 and b = 0.5 for p = 0, 1, 2, 3. 



p = 


WKB 


Time evolution 


L 


Re(fc) 


-Im(fc) 


Re(fe) 


— Im(fc) 





1.2889 


0.5506 


1.250 (3.0) 


0.4980 (9.6) 


1 


1.5047 


0.5876 


1.606 (6.7) 


0.4867 (17.2) 


2 


1.9638 


0.4812 


1.962 (0.092) 


0.4805 (0.15) 


p = l 


WKB 


Time evolution 


L 


Re(fc) 


-Im(fc) 


Re(fe) 


— Im(fc) 





1.0812 


0.4670 


1.042 (3.6) 


0.4498 (3.7) 


1 


1.3245 


0.4963 


1.604 (21.1) 


0.463 (6.7) 


2 


1.7264 


0.4301 


1.725 (0.079) 


0.4295 (0.13) 


p = 2 


WKB 


Time evolution 


L 


Re(fc) 


-Im(fc) 


Re(fe) 


— Im(fc) 





0.8714 


0.3911 


0.8346 (4.2) 


0.3926 (0.38) 


1 


1.1311 


0.4137 


1.161 (2.64) 


0.3803 (8.1) 


2 


1.488 


0.3754 


1.488 (0.013) 


0.3749 (0.13) 


p = 3 


WKB 


Time evolution 


L 


Re(fc) 


— Im(fc) 


Re(fe) 


— Im(fc) 





0.6633 


0.3202 


0.6376 (3.9) 


0.3279 (2.4) 


1 


0.9284 


0.3363 


0.9413 (1.4) 


0.3204 (4.7) 


2 


1.2489 


0.3162 


1.249 (0.0056) 


0.3157 (0.14) 



results near /3 — 1, as shown in Fig. |3j 

It is worth noticing that the frequency k show an almost scaling behaviour on functions 
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Table II: Fundamental quasinormal frequencies with a = 2 and b = 0.5 for p = 4,5, 6. 



p = 4 


WKB 


Time evolution 


L 


Re(fc) 


— Im(fc) 


Re(fe) 


— Im(fc) 





0.4632 


0.2514 


0.4449 (4.0) 


0.2555 (1.6) 


1 


0.7211 


0.2607 


0.7244 (0.46) 


0.2438 (6.4) 


2 


1.0081 


0.2512 


1.008 (0.012) 


0.2509 (0.13) 


p = 5 


WKB 


Time evolution 


I 


Re(fc) 


— Im(fc) 


Re(fe) 


— Im(fc) 





0.2825 


0.1828 


0.2697 (4.5) 


0.1990 (8.8) 


1 


0.5179 


0.1843 


0.5187 (0.16) 


0.1828 (0.83) 


2 


0.7690 


0.1804 


0.7691 (0.010) 


0.1802 (0.082) 


p = 6 


WKB 


Time evolution 


L 


Re(fc) 


— Im(fc) 


Re(fc) 


— Im(fc) 





0.3135 


0.05970 


0.1485 (52.6) 


0.1290 (116.1) 


1 


0.3608 


0.1154 


0.3616 (0.22) 


0.1150 (0.34) 


2 


0.5890 


0.1135 


0.5889 (0.021) 


0.1134 (0.042) 



of a -1 , as shown in Fig. HI That happens for the imaginary as well as for the real parts of 
k except for very small values of a. We found a different behaviour just in the case L = 2, 
n = 2, for the values of a < 2 near the extremal case a = b. No instability has been found. 
For higher dimensions the real and imaginary parts of the frequency decreases. An exception 
is the case L = 2, n = 2: the real part of the frequency increases in the range < p < 3 
and decreases for the others values of p, but the imaginary part decreases when p increases 
as for all others values of L and n that we considered in this work. We have found that for 
a given value of L increasing the overtone number n the frequencies become more damped, 
as we expected. 
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Table III: High overtone quasinormal frequencies with a = 2 and b = 0.5. 



p = 


P = 


= 1 


L 


n 


Re(fc) 


— Im(fc) 


Re(fc) 


— Im(fc) 


1 


1 


0.985828 


1.79911 


0.892835 


1.58205 


2 


1 


1.47092 


1.61706 


1.3581 


1.39048 


2 


2 


0.408755 


2.80627 


0.538849 


2.55582 


p = 2 


P = 


= 3 


Z, 


n 


Re(fc) 


— Im(fc) 


Re(fc) 


— Im(/c) 


1 


1 


0.798307 


1.34085 


0.693083 


1.09224 


2 


1 


1.22235 


1.18843 


1.0673 


0.990437 


2 


2 


0.638727 


2.20914 


0.690423 


1.82098 


p = 4 


P = 


= 5 


£ 


n 


Re(fc) 


— Im(fc) 


Re(fc) 


-Im(fc) 


1 


1 


0.572698 


0.841449 


0.439874 


0.587662 


2 


1 


0.895481 


0.781983 


0.710916 


0.55703 


2 


2 


0.681196 


1.40916 


0.609848 


0.980543 



p = 6 


L 


n 


Re(fc) 


— Im(fc) 


1 


1 


0.325148 


0.365934 


2 


1 


0.568875 


0.345275 


2 


2 


0.537292 


0.587244 



Although in general the calculation of the quasinormal frequencies can be only made using 
numerical methods, in the present scenario there is an important limit where an analytic 
expression is available. Expanding the effective potential in terms of small values of 1/L 
and using the WKB method in the lowest order (which is exact in this limit), we obtain: 

cu 2 = L 2 T(x m ) -i(n+^J LA(x m ) , (14) 
12 




Figure 3: Effect of (3 on the behavior of oj for p = 2 with a = 2, b = 0.5 and L = 0. Two 
different numerical methods were employed. They are consistent for small and large enough (3, but 
discrepant near (3 = 1. 



where T(x) = A(x) = -^y^ [In (A(x)/C(x))\ + T" . The peak of effective 

potential is determined by V(x)', and occurs at x m = [— 1c\j (c 2 + (c\ — Sci) 1 / 2 )] 1 ^ 7 p \ 
with ci = (7-p)(ab) 7 - p and c 2 = -(9 -p)a 7 - p . 

Far from the horizon the effective potential (with @ — 0), in terms of r*, assumes the 
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Figure 4: Effect of the a parameter on quasinormal frequency. The p-brane parameters of are 
b = 0.5, L = 0; and j3 = 0. 



form 



if < p < 6 



L(L + 1) + (2a - b) 
if p = 6 and L = 



In r* 

_4 



+ o 



In r< 



(15) 



In r* 

4 



L(L + 1) [^ + (2a-6)^] 

if p = 6 and L > 

With this effective potential, it is shown 
evolves, at late time, according to 



271 ] that an initial data with compact support 



(16) 



Therefore, at asymptotically late times the massless perturbation decay as a power-law tail. 

The power- law coefficient a(p,L) reflects the potencial asymptotic behavior. For p 
1, 3, 5, 6, a(p, L) can be analitycally determined using the results in 



2j: 



a(p, L) 



2L - p + 8 with p = 1, 3, 5 



(17) 



2L + 3 



with p = 6 



14 



For p = 0, 2, 4, our numerical resuts suggest a similar expression 

a(p,L) = 2L -p + 10 with p = 0,2,4 (18) 



The tails are confirmed by the time-dependent approach. We illustrate these results in Fig. 

El 




Figure 5: (Top) Tails for several values of p. The power-law coefficients estimated from the 
numerical data (with t > 250) are: -10.01 (p = 0), -5.07 (p = 3) and -3.15 (p = 6). The analytical 
results (indicated by straight lines) are: -10 (p = 0), -5 (p = 3) and -3 (p = 6). The p-brane 
parameters are a = 2, b = 0.5, L = and (3 = 0. (Bottom) Massive tail for p = 6. The envelope 
power-law coefficient estimated from the numerical data (with t > 9000) is —0.84. The analytical 
result (indicated by a straight line) is — 5/6 w —0.833. The p-brane parameters are a = 2, b = 0.5, 
L = and (3 = 1. 

In the massive case, the asymptotic form of the effective potential changes. For large r* 
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we have 



if < p < 5 



/3 2 +[/3 2 + L 2 + 2L + |]^ + 0(^ 
if p = 5 



(3 2 (l + ^)+o(4 
if p = 6 and L = 



(19) 



/3 2 (l + ^a) + [/3 2 6(& - a) 
+L(L + 1)]^ + 0(W) 

if p = 6 and L > 

We have observed from the numerical simulations that the late-time tail have the form 



^ L ~ sin(/3t)f 



-7(P.i) 



(20) 



If p = 6, the results in 28|-|30| apply, and the coefficient in the power-law envelope can be 
determined analytically: j(p = 6, L) = 5/6. This result is illustrated in Figj5j For other 
values of p the analytical problem remains open. 



V. EXTREME CASE 



The analysis of the extreme case geometry is more subtle. If p = 3, we have a black 
hole solution and the problem is clearly formulated. If p = 6, we have a naked time-like 
singularity and the Cauchy problem is not well-posed (without additional conditions at the 
singularity). This class of solution will not be treated in the present work. 

The novelty is the geometry with a null singularity. As discussed before, we have a well- 
posed initial value problem. We propose here to define the quasi normal modes in the same 
way they were defined in the black hole scenario. This definition will be justified considering 
the wave problem in the following. 
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The effective potential for the scalar field perturbation in the extreme case scenario is 
obtained by taking a = b in (j!2p . This potential looks similar to the non extreme case 
analog, and in terms of the tortoise coordinate, it tends to zero as r t -oo and — > oo, 
what implies that the effective one-dimensional wave problem is similar to the previous non- 
extreme case. A bounded perturbation will therefore decay in time, what justifies the quasi 
normal mode definition adopted. As a side remark, we observe that for p = 6 the potential 
diverges near the horizon, a consequence of the time-like nature of the singularity at r = a. 
We have computed the quasi normal frequencies for p < 5. The results are shown in Table 
IIVI We have sensible differences, by factors of order three. 

For L = 2, from n = to n = 1 we observe an increase in the decay rate. We found that 
the imaginary part increases, in the case p = from L — 2, n — 2 to n — 1, in contrast with 
the behaviour found in the non extreme case. Otherwise, results are very similar to the non 
extreme case. 



VI. FINAL REMARKS 



We studied the scalar perturbations of the full black p-brane solutions of ten dimensional 
type IIB Supergravity. The near horizon limit of extremal p— branes is an AdS( p+ 2) x S^ 8 ~ p ^ 
space-time, which is dual to a (p+l)-dimensional conformal field theory at zero temperature. 
If we have an event horizon, the near horizon limit is a (p + 2)-dimensional AdS black hole 
times a sphere S^ 8 ~ p \ dual to a field theory at finite temperture in (p + 1) dimensions. 
We obtained the same quasinormal spectrum using the standard procedure of considering 
a probe scalar field in the backgroud geometry with a gauge invariant formalism. The 
quasinormal modes structure in such a complex problem is amazingly simple. Allowing for 
a non vanishing separation constant, later related to the glueball mass, the result is also very 
simple, displaying an almost scaling behaviour. The tensor and scalar modes are exactly the 
same, leading to a simplicity of the results as well. Implications for the quark-gluon-plasma 
using the AdS/CFT relation awaits further analysis. 
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Table IV: Scalar quasinormal frequencies for the extreme case (a = b = 1). 



p = 


p=l 


L 


n 


Re(fc) 


-Im(fc) 


Re(fc) 


-Im(fc) 








2.49971 


0.955112 


2.04071 


0.847233 


1 





3.07066 


1.01771 


2.68902 


0.863596 


1 


1 


2.41319 


2.08326 


1.9794 


2.17509 


2 





3.88648 


0.931587 


3.47584 


0.802949 


2 


1 


3.15979 


2.7734 


2.88255 


2.40614 


2 


2 


0.0876808 


2.4072 


0.978905 


3.54121 


p = 2 


p = 3 


L 


n 


Re(fc) 


-Im(fc) 


Re(fc) 


-Im(fc) 








1.6804 


0.659786 


1.51662 


0.429531 


1 





2.35616 


0.685455 


2.09221 


0.522615 


1 


1 


1.79834 


1.90349 


1.7917 


1.46652 


2 





3.10546 


0.658724 


2.79837 


0.513388 


2 


1 


2.67848 


1.97719 


2.52534 


1.55801 


2 


2 


1.52044 


3.22954 


1.92728 


2.57206 



p = 4 


L 


n 


Re(fc) 


— Im(fc) 








1.40021 


0.32376 


1 





2.00824 


0.340953 


1 


1 


1.82987 


0.991697 


2 





13.399173 


15.362158 


2 


1 


2.52447 


1.05096 


2 


2 


2.19558 


1.73277 



18 



Acknowledgments 

This work has been supported by FAPESP and CNPq, Brazil as well as ICTP, Trieste. 



[1] J. Polchinski, String Theoy, vols. 1,2, Cambridge University Press, (1998); B. Zwiebach, A 

First Course in String Theory, Cambridge University Press, (2004). 
[2] J. Maldacena, Adv. Theor. Math. Phys. 2, 231, (1998). 

[3] O. Aharony, S. Gubser,J. Maldacena, H. Ooguri, Y. Oz, Phys. Reports 323, 183, (2000). 
[4] G. W. Gibbons and P. K. Townsend Phys. Rev. Lett. 71, 3754, (1993). 

[5] E. Shuryak, Progr. Part. Nucl. Phys. 53, 273, (2004); S. Gubser, Gen. Relat. Grav. 39, 1533, 

(2007); S. Gubser Phys. Rev. D74, 126005, (2006); J. Friess et al. JEEP 04, 080, (2007). 
[6] G. Policastro, D. T. Son, A. O. Starinets, JEEP 09, 043, (2002); D. T. Son, A. O. Starinets, 

JEEP 09, 042, (2002); D. T. Son, A. O. Starinets, Ann. Rev. Nucl. Part. Set., 57, 95, (2007); 

D. T. Son, A. O. Starinets, JEEP 03, 052, (2006). 
[7] P. Kotvun, D. T. Son, A. O. Starinets, Phys. Rev. Lett, 94, 111601,(2005); P. Kotvun, D. 

T. Son, A. O. Starinets JEEP, 0310, 064, (2003); A. O. Starinets, Phys. Rev., D66, 124013, 

(2002). 

[8] C. Csaki, H. Ooguri, Y. Oz and J. Terning, JHEP 9901, 017 (1999). 
[9] H. Boschi-Filho and N. R. F. Braga, JHEP 0305, 009 (2003). 
[10] A. S. Miranda, C. A. Ballon Bayona, H. Boschi-Filho and N. R. F. Braga, JHEP 0911, 119 
(2009). 

[11] S. Chandrasekhar, The Mathematical Theory of Black Eoles, New York, Oxford University 
Press, (1983). 

[12] H. Kodama, A. Ishibashi, Prog. Theor. Phys., 110, 701, (2003); H. Kodama, A. Ishibashi, O. 

Seto, Phys. Rev., D62, 064022, (2000). 
[13] K. D. Kokkotas, D. G. Schmidt, Liv. Rev. Relat, 2, (1999); E. W. Leaver, Phys. Rev. D45, 

4713, (1992); C. Gundlach, R. H. Price, J. Pullin, Phys. Rev., 49, 883, (1994); E. Abdalla, R. 

A. Konoplya, C. Molina, Phys. Rev., D72, 084006, (2005). 
[14] V. Cardoso, J. P. S. Lemos, Phys. Rev., D67, 084020, (2003); R. A. Konoplya, A. Zhidenko, 

JEEP, 0406, 037 (2004); R. A. Konoplya, Phys. Rev. D68, 124017, (2003). 



19 



[15] B. Wang, C. Y. Lin, E. Abdalla, Phys. Lett, B481, 79, (2000); B. Wang E. Abdalla, C. 
Molina, Phys. Rev., D63, 084001, (2001); C. Molina, D. Giugno, E. Abdalla, A. Saa, Phys. 
Rev., D69, 104013, (2004); E. Abdalla, B Cuadros Melgar, C. Molina, A. Pavan, Nucl. Phys., 
B752, 40, (2006). 

[16] Jaqueline Morgan, Vitor Cardoso, Alex S. Miranda, C. Molina, Vilson T. Zanchin, Journal of 

High Energy Physics, 0909, 117 (2009). 
[17] G. T. Horowitz, A. Strominger, Nucl. Phys. B360, 197, (1991). 
[18] A. Erdelyi, Higher Transcendental Functions, 2, 232, (1953). 
[19] Hans-Peter Nollert and Bernd G. Schmidt, Phys. Rev. D 45, p. 2617 (1992). 
[20] R. A. Konoplya, Phys. Rev. D 68, 024018 (2003). 
[21] C. Gundlach, R. Price, and J. Pullin, Phys. Rev. D49, 883 (1994). 

[22] P. R. Brady, C. M. Chambers, W. Krivan and P. Laguna, Phys. Rev. D 55, 7538 (1997). 
[23] C. Molina, D. Giugno, E. Abdalla, A. Saa, Phys. Rev. D69, 104013 (2004). 
[24] R. A. Konoplya and A. Zhidenko, Nucl. Phys. B 777 (2007) 182. 
[25] R. A. Konoplya and A. Zhidenko. larXiv: 0809 .20481 [hep-th], 
[26] R. H. Price, Phys. Rev. D5, 2419 (1974). 

[27] E. S. C. Ching, P. T. Leung, W. M. Suen and K. Young, Phys. Rev. D52, 2118 (1995). 

[28] H. Koyama, A. Tomimatsu, Phys. Rev. D 65, 084031 (2002). 

[29] H. Koyama, A. Tomimatsu, Phys. Rev. D 64, 044014 (2001). 

[30] R. A. Konoplya, A. Zhidenko and C. Molina, Phys. Rev. D 75, 084004 (2007). 



20 



